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"^ Abstract 

r^ We apply flicker-noise spectroscopy (FNS), a time series analysis method 

, Ph operating on structure functions and power spectum estimates, to study the 

clinical electroencephalogram (EEG) signals recorded in children/adolescents 
^ (11 to 14 years of age) with diagnosed schizophrenia-spectrum symptoms at 

■^ the National Center for Psychiatric Health (NCPH) of the Russian Academy 

of Medical Sciences. The EEG signals for these subjects were compared with 
the signals for a control sample of chronically depressed children/adolescents. 
1^ The purpose of the study is to look for diagnostic signs of subjects' suscep- 

O tibility to schizophrenia in the FNS parameters for speciflc electrodes and 

,__! cross-correlations between the signals simultaneously measured at different 

IJ points on the scalp. Our analysis of EEG signals from scalp-mounted elec- 

trodes at locations -F3 and F4, which are symmetrically positioned in the 
left and right frontal areas of cerebral cortex, respectively, demonstrates an 
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^ essential role of frequency-phase synchronization, a phenomenon represent- 
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ing specific correlations between the characteristic frequencies and phases of 
excitations in the brain. We introduce quantitative measures of frequency- 
phase synchronization and systematize the values of FNS parameters for the 
EEG data. The comparison of our results with the medical diagnoses for 84 
subjects performed at NCPH makes it possible to group the EEG signals into 
4 categories corresponding to different risk levels of subjects' susceptibility 
to schizophrenia. We suggest that the introduced quantitative characteris- 
tics and classification of cross-correlations may be used for the diagnosis of 
schizophrenia at the early stages of its development. 
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1. Introduction 

The objective diagnosis of psychiatric disorders, such as schizophrenia, 
at the early stages of their development is complicated by the lack of re- 
liable instrumental methods able to adequately describe these disorders at 
their onset [U El El II]- For example; the methods of electro- and magne- 
toencephalography are rarely used for the diagnosis of psychiatric disorders, 
though there are experimental studies suggesting a high potential of these 
methods O El [7]. This implies that clinical electroencephalogram (EEG) 
studies are of high concern for advancing the state of the art in the analysis, 
diagnosis, and prognosis of these disorders. 

Some aspects of the application of spectral methods to the analysis of 
EEG rhythmic activity in specific frequency ranges were previously discussed 
by Borisov et al. [8]. It was shown that one of the essential population signs 
of schizophrenia is the dampening of a-activity accompanied by the ampli- 
fication of low-frequency 6- and ^-activity. Segment analysis and statistical 
combinatorial analysis of spectral and segmental characteristics by pattern 
recognition were used to study the EEG temporal dynamics. The analysis 
demonstrated that the patients suffering from schizophrenia were character- 
ized by reduced amplitudes and lengths and increased variabilities of quasis- 
tationary a-activity segments, as compared to the controls. Kostyuchenko 
and Kaplan [9] observed in patients suffering from schizophrenia an increase 
in the spectral density of EEG fluctuations in the low-frequency range (be- 
low 0.25 Hz). A significant role of stochastic components, specifically flicker 



noise, in EEG signals was discussed in depth by West et al. [10], AUegrini 
et al. [n], Aquino et al. [12]. 

However, these conclusions are just empirical observations that cannot 
be clearly traced to the internal physiological mechanisms of psychiatric dis- 
orders. We believe there is an essential phenomenon of frequency-phase syn- 
chronization in the functional activity of some pathology-specific cortex areas 
that should be studied to provide a better understanding of the physiological 
aspects of schizophrenia as well as other neurodegenerative and psychiatric 
disorders. We expect that a cross-correlation analysis of frequency-phase 
synchronization will be more robust to the individual-specific variations of 
biomedical signals, which generally impede the conventional diagnosis. It 
should be noted that such an analysis cannot completely replace the charac- 
terization of specific EEG signals, but should rather complement it to develop 
a more comprehensive picture of the interactions between specific areas of 
the cortex. 

Frequency and phase synchronization, manifestation of specific correla- 
tions between characteristic frequencies and phases of the excitations in dif- 
ferent parts of the cortex (specific neural ensembles), and synchronization 
of the excitation amplitudes are the necessary conditions for the brain to 
function as an integral system [131 [IH [13 HB] • A normally functioning brain 
responds to external actions on the human organism by establishing some 
optimal level of such synchronization. A significant deviation from this op- 
timal level, such as an anomalously high level of synchronization or lack of 
synchronization, may be considered as an indicator of a pathology in brain 
activity. In the fields of medicine where digitized data, such as EEG [T71 [18] 
or magnetoencephalograms (MEG) [ini[2n], are already recorded for the di- 
agnosis, this optimal level can be estimated by analyzing the correlations 
between the signals recorded at spatially separated areas. 

Studies of cooperative neural processes in bioelectric brain activity, which 
bring about frequency-phase synchronization and other similar phenomena, 
are relevant not only for the diagnosis of psychiatric disorders, but also for 
some neurodegenerative disorders, such as Parkinson's disease [21] and tem- 
poral lobe epilepsy [22] , and are important for the understanding of human 
cognitive abilities [23]. The characteristic frequencies and relations between 
phase differences for EEG and MEG signals recorded at spatially separated 
electrodes are, as a rule, estimated using the Hilbert and Fourier transforms 
[T3| [21] or wavelet analysis [2S]. Recently, the methods of cross-correlation 
analysis were also introduced into the analysis of mutual dynamics of phys- 



iological signals [251 EZl ESI ES]- Panischev et al. PB] applied the formalism 
of memory functions to study stimulated neuromagetic responses, MEG sig- 
nals simultaneously recorded by different pairs of superconducting quantum 
interference device (SQUID) sensors located in spatially separated cortex 
areas. The analysis of synchronization effects in the mutual dynamics of 
neuromagnetic responses to flickering-color stimuli showed an increased level 
of frequency-phase synchronization in control (healthy) subjects compared 
to the patient suffering from photosensitive epilepsy. A detrended cross- 
correlation analysis (DCCA) method was proposed and applied to quan- 
tify cross-correlations in the presence of nonstationarity [2Z1 EO]- Another 
method, random matrix theory (RMT), which was developed to measure 
cross-correlations in collective modes for a large number of simultaneously 
recorded data [ST, '32], was employed to detect long-range power-law correla- 
tions in physiological data [28] . Bob et al. |29j reported that cross-correlations 
between pairs of EEG time series are inversely related to dissociative symp- 
toms (psychometric measures) in 58 patients with paranoid schizophrenia. 

In this study, we look for diagnostic signs of the susceptibility of children 
and adolescents to schizophrenia at its onset by analyzing the EEG signals 
recorded from scalp-mounted electrodes at locations F^ and F^. The clinical 
EEG signals were recorded in children/adolescents (11 to 14 years of age) 
with diagnosed schizophrenia-spectrum symptoms at the National Center 
for Psychiatric Health (NCPH) of the Russian Academy of Medical Sciences 
(Moscow, Russia). The EEG signals for these subjects were compared with 
the signals for a control sample of chronically depressed children/adolescents, 
also recorded at NCPH. 

We perform the analysis using flicker-noise spectroscopy (FNS), a time 
series analysis method that introduces information parameters characterizing 
the components of stochastic signals in different frequency ranges and two- 
parameter cross-correlation functions [331 EH ESI ESI EZl EH EH] • The FNS 
method separates the analyzed signal into three components: low-frequency 
regular component corresponding to system-specific "resonances" and their 
interferential contributions, stochastic random-walk component at larger fre- 
quencies corresponding to dissipation effects, and highest-frequency inertial 
"spike" component corresponding to flicker noise [iQl chap. 4.3]. It should 
be noted that the latter is taken into consideration due to the intermittent 
dynamics on various space-time scales attributed to complex (multiparti- 
cle, nonlinear) interactions, dissipation processes, and inertia characterizing 
thermodynamically open complex systems [H]. 



The FNS method was previously used to study frequency-phase synchro- 
nization between spatially separated cortex areas for the case of photosen- 
sitive epilepsy [371 EH] . Neuromagnetic responses to external flickering-color 
stimuli, the prolonged action of which can cause epileptic seizures, were ex- 
amined. The FNS analysis detected a disorder-specific feature, manifestation 
of high-frequency components (approximately 50 and 100 Hz) in the power 
spectum estimates for some SQUID-sensors. The analysis of FNS cross- 
correlation functions showed substantial disruptions in frequency-phase syn- 
chronization for the MEG dynamics of patient cortex areas compared to the 
synchronization level of control subjects [S^ . 

The paper is structured as follows. In Section 2, we provide the fundamen- 
tals of FNS and present the parameterization algorithm and two-parameter 
cross-correlation function. In Section 3, we perform the parameterization of 
EEG signals at electrodes F3 and F4. In Section 4, we present the results 
of cross-correlation analysis of the signals, which suggest that the level of 
frequency-phase synchronization can be considered as a diagnostic sign of 
subjects' susceptibility to schizophrenia. In Section 5, we introduce quanti- 
tative measures of frequency-phase synchronization that allow us to partition 
all 84 subjects of the clinical study into 4 categories corresponding to differ- 
ent levels of the risk for developing schizophrenia. In Section 6, we discuss 
the potential capabilities of frequency-phase synchronization analysis for the 
diagnosis of neurological and psychiatric disorders. 

2. Principles of flicker-noise spectroscopy and basic relations 

Here, we will only deal with the basic FNS relations needed to under- 
stand the parameterization procedure and cross-correlation function. FNS is 
described in more detail elsewhere [331 EH ESI ESI EH] • 

In FNS, all introduced parameters for signal V{t), where t is time, are 
related to the autocorrelation function 

^(r) = (r(t)V(t + r))^^,, (1) 

where r is the time lag parameter (0 < r < Tm) and Tm is the upper bound 
for r {Tm < T/2). This function characterizes the correlation in values of 
dynamic variable V at higher, t + r, and lower, t, values of the argument. 
The angular brackets in relation ([I]) stand for the averaging over time interval 
T-t: 

{{..■))r-r = ^ r \-)dt. (2) 



The averaging over interval T — t implies that all the characteristics that 
can be extracted by analyzing functions iP{t) should be regarded as the 
average values on this interval. To extract the information contained in ^(r) 
{{V{t)) = is assumed), the following transforms, or "projections", of this 
function are analyzed: cosine transforms (power spectrum estimates) S{f), 
where / is the frequency, 

S{f) = 2 / {V{t)V{t + ti))T-. cos(27r/ti) dt, (3) 

Jo 

and its difference moments (Kolmogorov transient structure functions) of the 
second order $*^^)(r) 

^^'\T) = {[V{t)-V{t + T)f)^^^. (4) 

The information contents of S{f) and $'^^^(r) are generally different, and 
the parameters for both functions are needed to solve parameterization prob- 
lems. By considering the intermittent character of signals under study, in- 
terpolation expressions for the stochastic components $s^^^(r) and Ss{f) of 
S{f) and $'^^^(r), respectively, were derived using the theory of generahzed 
functions by Timashev [33] . It was shown that the stochastic components of 
structure functions $*^^)(r) are formed only by jump-like irregularities ("ran- 
dom walks"), and stochastic components of functions S{f), which character- 
ize the "energy side" of the process, are formed by spike-like (inertial) and 
jump-like irregularities. 

2.1. Signal parameterization 

In FNS parameterization, the original signal V{t) is separated into three 
components: system-specific "resonances" and their interferential contribu- 
tions at lower frequencies, stochastic jump-like ("random walk") component 
at larger frequencies, and stochastic spike-like (inertial) component in the 
highest frequency range. For simplicity, we will further refer to jump-like 
and spike-like irregularities as "jumps" and "spikes", respectively. 

Let us write the basic interpolation expressions for stochastic components. 
The parameters characterizing the dynamic correlations on every level of the 
evolution hierarchy are assumed to be the same. Consider the simplest case, 
in which there is only one characteristic scale in the sequences of spikes and 
jumps ^[31]: 

$f (r) ^ 2a' [1 - r -\H^) ■ T (H,, r/T,)] ' , (5) 



where r(s,x) = Jexp(— t) ■ f^ ^dt and T{s) = r(s,0) are the complete and 

X 

incomplete gamma functions, respectively (x > and s > 0); a is the stan- 
dard deviation of the measured dynamic variable with dimension [V^]; Hi is 
the Hurst constant, which describes the rate at which the dynamic variable 
"forgets" its values on the time intervals that are less than the correlation 
time Ti. 

For asymptotic cases, we obtain the formulas: 

<l>f\T) = 2T-\l + Hi)-a'(^] , ^«1; (6) 



$(2)f^) = 2a2 



, Hi-l 
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l-T-^m-[^] expl-^^ 
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^ » 1- (7) 



The interpolating function for power spectrum component Sgsif) formed 
by spikes can be written as: 

Here, Sss{0) is the parameter characterizing the low-frequency limit of Sssif) 
and uq describes the degree of correlation loss in the sequence of spikes on 
the time interval Tq. 

The interpolating function for the power spectrum component Scj{f) 
formed by jumps is written as: 



SsjjO) 
(2vr/ri 



SsAf) - . , .::?:' ■ (9) 



where Ssj{0) is the parameter characterizing the low-frequency limit of Ssj{f)- 
Although the contributions to the overall power spectrum Ss{f) given by 
Eqs. (Is]) and (|9| are similar, the parameters in these equations can be much 
different: Ssj{0) 7^ Sss{0), Ti 7^ Tq, and 2Hi + 1 y^ Hq. This implies that the 
parameters in the expressions for the power spectrum and structure function 
generally have different information contents when the experimental time 
series V{t) is analyzed. In cases when the contributions of spikes and jumps 
into the overall stochastic component are comparable, which corresponds to 



relatively small differences in the values of parameters Uq and 2Hi + 1, it is 
convenient to use the following interpolating function for Ss{f)'- 

where Ss{0) and n are phenomenological parameters. 

The dynamics of complex systems includes both stochastic components, 
i.e., spikes and jumps, and system-specific slowly varying regular components 
associated with a set of frequencies. These frequencies correspond to internal 
and external resonances and their interferences. It should be noted that 
the whole set of resonance and interferential frequencies may get rearranged 
during the evolution of an open system. All the specific frequencies and their 
interferential contributions, which manifest themselves as oscillations in the 
dynamic variable V{t), will be be further called "resonant". 

It is assumed that signal V{t) can be presented as a linear superposition 
of stochastic component Vs{t) and resonant component K(t): 

Vit) = Vsit) + Vrit). (11) 

In this case, the autocorrelation function and power spectrum can be 
approximated as [311 EEl [37] : 

^{T)=^,{T) + Mr), (12) 

S{f) = Ss{f) + Sr{f). (13) 

Here, we assume that the resonant components are statistically stationary 
(they depend only on time lag r). This allows us to estimate ^ipriT) as an 
"incomplete" cosine transform of Sr{f) by applying the Wiener- Khinchin 
Theorem: 

/max 

Mr) ^ 2 j Srif) cos{2n fT)df, (14) 



where /max = 0.5/^, fd is the sampling frequency. It should be noted that 
Eq. ( 14 ) is an approximation applied to a finite discrete time series assuming 
the wide-sense stationarity of the resonant signal component. The resonant 
component $r (t) in this case is found by 

^^^\T) = 2[M0)-Mr)]- (15) 



The stochastic component of $'-^)(r) can then be estimated as 

$f(r)^$(2)(r)-$(2)(r). 



(16) 



Equations (12)-(16) allow one to sequentally separate out resonant and 



stochastic components of structure functions and power spectrum estimates 
for experimental time series and perform the parameterization of the com- 
ponents, which is demonstrated in Refs. [SU |38] . The complete FNS param- 
eterization algorithm in discrete form is listed in Appendix A. 

In total, six stochastic FNS parameters are introduced (Tqi, Ss(Tq^^), n, a, 
Ti, Hi): Tqi, the correlation time for jump- and spike-like irregularities after 
which the self-similarity observed in power spectrum estimate breaks down; 
Ss(Tq^^), the "spikiness" factor - power spectrum estimate at frequency Tq"\ 
which accounts for the "intensity" of jump- and spike-like irregularities in 
the highest-frequency interval; n, the flicker-noise parameter, which charac- 
terizes the rate of loss of correlations in the series of high-frequency irregu- 
larities in time intervals Tqi; a, the standard deviation of the value of the 
measured dynamic variable from the slowly varying resonant (regular) com- 
ponent, which is based solely on jump-like irregularities; Ti, the correlation 
time for jump-like irregularities in stochastically varying signal V{t)] Hi, the 
Hurst exponent (this estimate of the Hurst component is also referred to in 
literature as the Hausdorff exponent), which describes the rate at which the 
dynamic variable "forgets" its values on the time intervals that are less than 
Ti. 

2.2. Cross- correlation function 

The information about the dynamics of correlations in variables Vi{t) and 
Vj{t), measured at different points i and j, can be extracted by analyzing 
the temporal variations of various cross-correlation functions. Here, we will 
use the simplest "two-point" correlation expression characterizing the links 
between Vi(t) and Vj(t) 



qij{r,9i 



V,it)-V,it + T) 



.(2) 



^r' (r) 



V,it + e,j)-Vjit + 9ij + T) 



.(2) 



^r (r) 



T-T-\9ij\ 

(17) 
where r is the "lag" time (r > 0), 6ij is the "time shift" parameter. 

The cross-correlation expression qij{T]6ij) is a function of temporal pa- 
rameters r and 6ij, which can be represented as a three-dimensional plot. Of 
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most interest for the analysis are the intervals of r and 9ij where the cross- 
correlation function qij approaches positive unity (maximum level of positive 
correlations) or negative unity (maximum level of negative correlations). The 
value of Oij corresponding to maximum values of cross- correlation qij{T]dij) 
characterizes the cause-and-effect relation ( "flow direction" ) between signals 
Vi{t) and Vj{t). When 6^ > 0, the flow moves from point i to point j, when 
9ij < 0, from j to i. When the distance between points i and j is fixed, the 
value of 6ij can be used to estimate the rate of information transfer between 
these two points. 



The magnitude and behavior of the two-parameter expression (17) may 



significantly depend on the value of selected averaging interval T and upper- 
bound values of r and 9ij, which we will refer to as Tmax and 6'max- From the 
statistical reliability point of view, we set a constraint of rmax + |^max| < T/2. 
For conciseness, from now on we will refer to 6ij as 6. 

3. Parameterization of EEG signals at F^ and F4 electrodes 

The EEG signals were recorded using Ag/AgCl electrodes (BrainAmp, 
Brain Products GbmH, Germany) in wakeful relaxed children/adolescents 
with the eyes closed at 16 electrode sites (Oi, O2, -P3, -P4, Pz, T5, Tg, C3, 
C4, Cz, T3, T4, F3, F4, Fj, and Fg) set according to the international 10-20 
system. The participants were seated comfortably in a dimly lit, electrically 
shielded room. The electrodes were monopolarly referenced to coupled ear 
electrodes. Electrode impedances were monitored and kept below 10 KQ 
during the entire experiment. The EEG signals were continuously recorded 
at a sampling rate of 256 Hz and then down-sampled off-line to 128 Hz (/^ 
= 128 Hz). Only artifact-free EEG segments were used for the analysis. For 
every electrode, 7680 values of the electric potential were examined. 

Several longitudinal studies of schizophrenia suggest the existence of a 
cerebral degenerative process in, at least, some patients |l2l HSl HI] . Accord- 
ing to studies |l3l ESI HB] , this neurodegeneration may be most pronounced in 
the frontal lobes of the cerebrum. For this reason, we chose frontal electrodes 
F3 and F4 for our analysis. The examined signals Vp^it) and VpAit) for elec- 
trodes F3 and F4, respectively, were recorded in 39 healthy subjects (marked 
with prefix "S") and 45 subjects susceptible to schizophrenia. This initial 
grouping was done based on the medical diagnosis by specialists, which was 
performed at NCPH. Our analysis included the determination of subject- 
specific FNS parameters for the F3 and F4 signals and examination of cross- 
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correlations between the signals by Eq. (17). This comparative study allowed 
us to partition the EEG signals of 84 subjects into 4 categories corresponding 
to different levels of risk for developing schizophrenia, starting with group I 
(lowest risk) and ending with group 4 (highest risk). The subjects belong- 
ing to the first two groups may be characterized according to the medical 
diagnosis as "healthy" while the last two groups correspond to "susceptible 
to schizophrenia". For conciseness, in this and next sections we discuss in 
detail only selected examples of the analysis to illustrate our partitioning 
logic. The complete classification of all subjects is presented in Section 5. 

First, we would like to note that the FNS parameterization alone of EEG 
signals generated in the cortex areas under electrodes F^ and F4 was insuf- 
ficient to partition all 84 subjects into risk groups due to a high individual- 
specific variability of FNS parameters for the stochastic components of the 
signals. To unambiguously partition all subjects into 4 risk groups, we had 
to complement the FNS parameterization of F3 and F4 signals with their 



cross-correlation analysis using expression (17). 

The degree of scatter in the values of FNS parameters can be seen from 
Table [T| which shows the values of all FNS parameters for typical subjects 
of each of the four groups. The FNS parameterization was run using the al- 
gorithm listed in Appendix A for the averaging interval T equal to the series 
length (7680 points). Some typical dependences used in FNS parameteriza- 
tion are illustrated in Figs. [l]j4j The complete partitioning of all subjects into 
four groups corresponding to different risk levels for developing schizophrena 
together with the values of group-determining FNS parameters and quantita- 
tive estimates of cross-correlations is presented in Table[2j Table [T] shows that 
the values of FNS parameters for subjects assigned to different groups may 
significantly overlap, which implies that FNS parameterization alone cannot 
be used for the classification. The combined analysis of FNS parameters and 
cross-correlations discussed in Sections 4 and 5 allowed us to suggest the 
spikiness factor Ss{Tq^^) as the most functionally significant FNS parameter 
for the classification. Other FNS parameters listed in Table [T] have a rather 
limited value for the partitioning: it is possible to draw conclusions only on 
the form of their variability for subjects of different categories. 

The first two groups are denoted as "healthy" (group I) and "condition- 
ally healthy" (group II). It can be seen that the EEG signals of groups I 
and II are characterized by a relatively rapid loss of correlations in stochastic 
components, which follows from high values of the Hurst component {Hi > 1) 
and flicker-noise parameter n, and small correlation times Ti and Tqi (Figs. 
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Table 1: FNS parameters for typical members of each group, 
strong nonstationarity in EEG signals 



(*) denotes the cases of 



Subject 


Elcctr. 


O", 


Hi 


^1' fl' 




n 


T-oi, f^' 


Group I 


S47 


F:i 


360 


1.22 


2.30 


3.8 X lO'^ 


3.36 


2.64 


Fa 


324 


1.04 


2.45 


5.7 X lO'^ 


2.98 


2.58 


S163 


F-A 


145 


1.87 


0.83 


42.8 


4.45 


1.30 


Fa 


158 


2.03 


0.86 


61.8 


4.43 


1.48 


S165 


F-i 


293 


1.84 


1.46 


801.26 


3.96 


2.39 


Fa 


263 


1.68 


1.44 


762.35 


3.81 


2.19 


Group II 


S177* 


Fi 


209 


0.02 


3.2 X 10'' 


157 


4.10 


1.30 


Fa 


170 


0.05 


9.1 X lO'' 


24 


4.35 


0.87 


156 


F; 


282 


1.13 


4.64 


2.42 X 10" 


2.58 


6.85 


Fa 


256 


1.01 


5.61 


2.7 X 10* 


2.46 


7.68 


S42 


Fa 


196 


2.24 


0.57 


26.27 


4.95 


1.04 


Fa 


151 


2.77 


0.37 


12.68 


4.98 


0.86 


Group III 


575 


Fa 


291 


1.10 


2.60 


6.2 X lO'^ 


2.90 


3.00 


Fa 


235 


1.31 


1.86 


2.4 X lO'^ 


3.11 


2.51 


S31 


Fa 


382 


1.12 


3.24 


1.5 X 10** 


2.89 


3.95 


Fa 


415 


0.85 


5.45 


2.7 X 10** 


2.73 


4.96 


S12 


Fa 


133 


1.95 


0.75 


36.00 


4.43 


1.23 


Fa 


118 


3.02 


0.38 


32.55 


4.45 


1.18 


Group IV 


221 


Fa 


558 


0.56 


35.9 


6.04 X 10-' 


2.20 


26.5 


Fa 


520 


0.56 


37.9 


6.58 X 10=^ 


2.10 


29.8 


387_03 


Fa 


356 


1.28 


3.65 


1.2 X 10" 


3.07 


4.89 


Fa 


311 


1.35 


3.21 


7.8 X lO'^ 


3.11 


4.56 


573 


Fa 


232 


0.33 


37.15 


1.7 X lO'' 


1.72 


22.7 


Fa 


228 


0.46 


12.25 


5.9 X 10** 


1.87 


9.51 
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|3^, |3|d). The group of "conditionally healthy" subjects also includes 5 sub- 
jects (S59, S170, S177, 351, S43) with a strong nonstationarity in the EEG 
signals: the steady-state variance is not reached within the analyzed inter- 
val (formally, this corresponds to i^i <^ 1 and Ti :§> T). The subjects in 
groups I and II are also characterized by moderate values of spikiness fac- 
tor Ss(Tq^^) < 3 X lO'^fiV'^fJ^. The third and fourth groups are composed of 
subjects "susceptible to schizophrenia" and "highly susceptible to schizophre- 
nia" , respectively. The EEG signals for these categories are characterized by 
a high degree of correlations in the stochastic irregularities of both types, 
which can be seen from the lower values of Hi and n as compared to groups 
I and II, and high values of spikiness factor Ss(Tq^^) (up to 7 x lO^^V'^f^^). 
The correlation times Ti and Tqi are small (Figs, pp, pH) in this case as well. 
A significant scatter in the values of characteristic parameters was also 
observed for low-frequency spectra of the signals, i.e., in the peaks corre- 
sponding to a series of "resonances". These "resonance" frequencies can be 
attributed to individual-specific resonances and interferential effects caused 
by the interaction of intrinsic and external frequencies. The power spectrum 
estimates of EEG signals recorded by electrodes F^ and F^ for the subjects 
of group I (Fig. flk) display clear peaks in the frequency range from 5 to 12 
Hz (Fig. [2k). These peaks can be attributed to the brain a-rythm, which 
manifests itself at the times of rest, relaxation, or slight meditation with the 
eyes closed. The dynamics of the difference moment for these EEG signals 
has a small-scale structure (Fig. Ilk). The power spectrum estimates of EEG 
signals for the subjects of group II (Fig. ^p) contain more noise on the back- 
ground of which there are peaks corresponding to brain rythms (Fig. [2b). As 
compared to group I, the magnitude of the peaks is reduced and the peaks 
themselves are shifted to lower frequencies. These results are in agreement 
with prior results reported by Borisov et al. |8], where the dampening of 
a-activity and amplification of low-frequency S- and ^-activity is considered 
as a diagnostic sign of schizophrenia. The scale of fluctuations in the dif- 
ference moment is increased compared to the group I case (Fig. [4Jd). The 
EEG signals for the group III subjects (Fig. [Ik) are characterized by further 
reduction of the activity in the range from 5 to 12 Hz and higher peaks in the 
range of ultralow frequencies (0.4 — 4 Hz). The scale of fluctuations in the 
difference moment gets further increased (Fig. [Ik). In the group IV subjects 
(Fig. liji), the ultralow frequency dynamics becomes even more pronounced: 
peaks in the frequency range higher than 2.5 Hz have no noticeable effect and 
the large-magnitude peaks in the ultralow frequency range play the dominant 
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Figure 1: (Color online) EEG signals at electrode F3 for typical members of groups I (a), 
II (b), III (c), and IV (d). 
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Figure 2: (Color online) Linear-scale power spectrum estimates of EEG signals in Fig. [T] 
for typical members of groups I (a), II (b), III (c), and IV (d). 
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Figure 3: (Color online) Log-log-scale power spectrum estimates of EEG signals in Fig. IT] 
for typical members of groups I (a), II (b), III (c), and IV (d); 1 (blue) - power spectrum 
estimate for the experimental data, 2 (red) - stochastic interpolation Ss{f). 
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Figure 4: (Color online) Structure functions of EEG signals in Fig. [Tlfor typical members 
of groups I (a), II (b), III (c), and IV (d); 1 (blue) - structure function for the experimental 
data, 2 (red) - resonant component of the structure function, "cross" (green) - total FNS 
interpolation. 
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role (Fig. |2]i). This observation is in agreement with the results of study [H], 
where it was noticed that the subjects suffering from schizophrenia have an 
increased spectral density of EEG signal fluctuations in the low-frequency 
range. In the group IV subjects, the structure of the difference moment is 
dominated by large-scale fluctuations (Fig. Ill: 



4. Frequency-phase synchronization in EEG signals at F^ and F4 
electrodes 

We analyzed the dynamics of cross-correlations between EEG signals at 



electrodes F3 and F4 using Eq. (17). The two-parameter cross-correlation 
dependences were calculated for adjacent intervals with length T = 800/^^ = 
6.25 s for all 84 subjects. This implies that the studied time interval 7200 fj^ 
= 56.25 s was partitioned into 9 subintervals (1...800) /^"^ ,..., (6401. ..7200) 
f^^, for each of which two- and three-dimensional plots of the cross-correlation 
function were built. The ranges of the input parameters were < r < r^ax = 
O.IT and — 6'max < 6* < 6'jnax5 where ^max = 0.4T. Here, we present plots only 
for selected subintervals to illustrate the general structure of observed cross- 
correlation variations. 

Figures [5]J6] show three-dimensional plots of cross-correlation function 
g(r, 6) for the F3 and F4 EEG signals of some of the typical subjects that 
were considered in Section 3. These examples support our working hypothe- 
sis stating that the degree of disruption in frequency-phase synchronization 
within specific cortex areas (in this case, left and right frontal areas) can be 
used as a diagnostic sign of subjects' susceptibility to schizophrenia-spectrum 
disorders. 

Figure |5] illustrates the cross-correlation function and its plane projec- 
tions for the EEG signals of subjects S47 (Figs. [5^, ^) and S177 (Figs. 
|5p, |5Ji) built for one of the subintervals of length 800 f^^. Frequency- 
phase synchronization between the signals F3 and F4 recorded in subject 
S47 manifests itself as a periodic dependence of g(r, 6) on 6, which can 
be best illustrated by looking at the cross-section of the three-dimensional 
plot at Tq = const. Figure pp shows the cross-section at Tq = 40/^^ (this 
choice of cross-section will be explained later). The maximum value of 
cross-correlation function g(40/^\0) ~ 0.89 suggests a high level of cor- 
relation between simultaneously recorded Vp^ and Vfa- A local maximum 
g(40/^^, —15/^7^) ~ 0.35 points to a correlation between Vp^ and Vf4, shifted 
by 9i = 15/^^ backwards. 9i = 15/^^ corresponds to a phase shift of 



Figure 5: Typical 3-D plots of cross-correlation function q{T, 9) for EEG signals at 
locations F3 and F/^ and their cross-sections at tq = 40/^^ built for subject S47, group I 
(a, c) and S177, group II (b, d). 
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Figure 6: Typical 3-D plots of cross-correlation function q{T, 9) for EEG signals at 
locations F3 and i^4 and their cross-sections at tq = 40/^^ built for subject 575, group III 
(a, c) and 221, group IV (b, d). 
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Stt/i^i = 27r, where /i = 1/15/^ ~ 0.5 Hz. This imphes that signal VF4(t) 
"lags" from signal Vps^t) by phase 2tt. Likewise, if we consider a local maxi- 
mum g(40/^^, 15/^^) ~ 0.25, we can conclude that signal VFi{t) goes ahead 
of signal VF3{t) by the same time interval or phase 27r. Therefore, we ob- 
serve a phase synchronization of signals Vps and Vp^, which takes place in 
the cortex of subject S47 within studied time interval. The presence of local 
maxima g(40/^^, —30/^^) and g(40/^^, 30/^^) in Fig. |5p with values exceed- 
ing 0.1 implies that the phase synchronization of signals Vps and Vfa occurs 
with a constant frequency /i, which allows us to use the term "frequency- 
phase synchronization". It should be noted that the observed frequency /i 
is individual-specific and can vary near its base value. 

The threshold value of 0.1, which was marked on the plane projections 
of cross-correlation function q{TQ, 6) by a horizontal dashed line, was selected 
based on the condition of closest agreement between our cross-correlation 
analysis and medical diagnosis at NCPH. The selection of tq = 40/^^ for 
the cross-section of cross-correlation function q{TQ,9) was done based on the 
observed characteristic period ~ 15/^^ in recorded signals. It can be assumed 
that at time lags tq exceeding two periods of oscillations the correlation links 
in stochastic components of difference moments are lost and the steady-state 
variance is reached [311 ES] . 

The cross-correlation dependences qij.O) and g(4O/^^,0) for the F3 and 
F4 signals of subject SI 77 also demonstrate the frequency-phase synchroniza- 
tion. But in contrast to S47, here, the number of local maxima exceeding 
0.1 is reduced and their shape is significanly deformed. In other words, the 
effects of frequency-phase synchronization are less pronounced in the case of 
group II as compared to group I. At the same time, a comparative analysis 
of frequency-phase synchronization and FNS parameters with medical diag- 
nosis still indicates that the susceptibility to schizophrenia for group II is 
relatively low. 

A similar logic applies to the analysis of cross-correlation dependences 
g(r, 6) and their plane projections g(40/^^, 0) of the EEG signals for groups 
III and IV (575, 221) (Fig. |6]). In general, the structure of three-dimensional 
plots for cross-correlation function gets largely deformed, and the number 
and magnitude of local maxima substantially decrease for groups with higher 
identification numbers (III and IV). 
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5. Classification of EEG signals based on degree of frequency-phase 
synchronization 

The results presented above suggest that quantitative estimates charac- 
terizing the degree of frequency- phase synchronization in signals Vp^lt) and 
VF4:it) should be introduced in addition to FNS parameters for the signals 
themselves in order to objectively and unambiguously estimate subjects' sus- 
ceptibility to schizophrenia. To estimate the degree of synchronization, we 
can use the plane projection of g(r, 9) at value tq and count the number of 
pairs of local maxima along the 6 axis, which are symmetric with respect to 
^ = 0, in the selected range \9\ < 150 f^^ for each sub interval with length 
T = 800/^^ = 6.25 s. The selected ranges \9\ < 150/^^ with respect to 
^ = are marked by vertical dashed lines in Figs. [5]j6| The magnitudes 
of all counted local maxima must exceed a specific threshold value (0.1 in 
this case). The number of such pairs n^ per unit time averaged for all 9 
sub intervals can be used to characterize the frequency fs = Us/T of "deep 
synchronization" for signals Vp^it) and Vf4(^) and be regarded as the basic 
criterion in assigning a subject to a specific category. The discreteness of 
pair number n^ implies that the values of fs are multiples of 0.16 Hz for the 
chosen value of T. Therefore, the determined values of frequency fs should 
only be regarded as rough estimates ("rounded" up to 0.16 Hz). At the same 
time, the number of observed pairs Hs can be considered as an objective crite- 
rion for determining the susceptibility of a child/adolescent to schizophrenia 
(lower risk for developing schizophrenia in adult age) for the given value of 
T. 

Our analysis showed that for 23 subjects out of total 84 the frequency fs 
was at least 0.48 Hz. For 2 subjects (S47 and S153) this frequency was 1.28 
Hz, and for subject S173 it reached the value of 1.6 Hz. These 3 values signif- 
icantly exceed the average fs for the group of subjects labeled as "healthy" 
by NCPH. These anomalies in the level of frequency-phase synchronization 
require an additional analysis and may be related to possible precursors of a 
hidden disorder. There are known examples in medicine when the increase 
of synchronization level in the dynamics of specific cortex areas beyond some 
threshold value is considered as an indicator of a disorder. An example of an 
anomalously high activity in neural ensembles associated with an epileptic 
seizure is discussed in Refs. [2S1 EH] • It is possible that a similar event took 
place in subject S47 on some of the considered time subintervals. As the 
authors do not have any information on the later examinations of these 3 



22 



Table 2: Classification of subjects into categories corresponding to different risk levels for 
developing schizophrenia based on the analysis of F3 and F4 EEG signals: I - "healthy" , 
II - "conditionally healthy" , III - "susceptible to schizophrenia" , IV - "highly susceptible 
to schizophrenia" . (*) denotes the cases of strong nonstationarity in EEG signals 



I 


II 


III 


IV 


Subject 


/„Hz 


Xl03, 

/^vV/d 


Subject 


/s,Hz 


S.(Tor) 
xIQS, 

MVV/d 


Subject 


/s,Hz 


xlO'^, 
MVV/d 


Subject 


/s,Hz 


xl03, 

/^vV/d 


S94 


0.48 


0.26 


416 


0.32 


0.01 


S12 


0.16 


0.04 


113 


<0.01 


575 


454 


0.48 


0.38 


S42 


0.32 


0.26 


586 


0.16 


0.08 


221 


<0.01 


658 


S182 


0.48 


0.68 


425 


0.32 


0.28 


S158 


0.16 


0.1 


307 


0.16 


7.2 


S165 


0.48 


0.8 


S154 


0.32 


0.54 


84 


0.16 


0.13 


S18 


0.16 


7.2 


S20 


0.48 


1.67 


S50 


0.32 


0.55 


83 


0.16 


0.13 


719 


0.16 
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S60 


0.48 
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0.32 


0.6 


S72 


0.16 


0.16 


387_03 


0.16 


11.6 
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0.64 


0.06 


S152 


0.32 


0.8 
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0.16 
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314 
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subjects labeled as "healthy" by NCPH, we will assign them to category I 
("healthy") in the classification described below. For 23 subjects, the fs pa- 
rameter was equal to 0.32 Hz; for 36, fs = 0.16 Hz; and for the two remaining 
subjects, fs < 0.01 Hz. The lower values of fs generally corresponded to a 
higher risk for developing schizophrenia (determined by medical diagnosis at 
NCPH). 

A comparative analysis of the values of Ss(Tq^^) (maximum value for 
signals F^ and F4) for the EEG signals belonging to different categories 
suggested that the spikiness factor can be used together with fs as a ba- 
sic diagnostic sign for schizophrenia-spectrum disorders. The higher values 
of 5*5 (Tq^^) and lower values of fs correspond to higher risks for develop- 
ing schizophrenia. As a result, the classification presented in Table [2] was 
built using the optimality criterion for the values of these two parameters 
that leads to the minimal discrepancy between the results of medical diag- 
nosis and proposed partitioning. The limiting cases of the partitioning are: 
(1) large values of parameter fs accompanied by small values of parameter 
Ss(Tq^^) correspond to a healthy subject; (2) small values of fs with large 
values of Ss(Tq^^) correspond to a high risk for developing schizophrenia. In 
contradictory cases, such as S26, S43, and 156, other FNS parameters were 
also used to assign the subjects to appropriate groups. For example, the high 
levels of "deep synchronization" {fs ~ 0.5) and stochasticity in subjects S26, 
S43, and 156 allowed us to conclude that the subjects belong to group II. 

The complete algorithm for partitioning subjects into 4 groups corre- 
sponding to different risk levels for developing schizophrenia is written as: 

(1) Group I. Necessary condition: /, > 0.48 Hz. If /, = 0.48 Hz and Ti < T, 
the following condition should also hold: Ss{Tq^^) < 3 x 10^/iV^/^^ If 
fs = 0.48 Hz and Ti > T (nonstationary), the following condition should 
also hold: ^.(Tq-/) < 6 x lO^VVd"^- If fs = 0.64 Hz, the following 
condition should also hold: ^^(^01^) < 3 x lO'^fxV^f^\ 

(2) Group II. Necessary condition: 0.32 < fs < 0.64 Hz. If /, = 0.32 
Hz and Ti < T, the following condition should also hold: Ss{Tqi^) < 
6x lOVVVd"^- If fs = 0.32 Hz and Ti > T (nonstationary), the following 
condition should also be true: 5,(Toi^) < 3 x lO^VVd"^- H fs = 0-48 Hz 
and Ti < T, the following condition should also hold: 3 x lO^^V'^f^^ < 
5,(Toi^) < 5 X lOVVVd"^- If fs = 0.48 Hz and Ti > T (nonstationary), 
the following condition should also hold: 6 x 10^/iV^/^^ < Ss{Tq^^) < 
5 X lO^'/iV^/^^. If fs = 0.64 Hz, the following condition should also hold: 
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(3) Group III. Necessary condition: 0.16 < fs < 0.32 Hz. If /, = 0.16 
Hz and Ti < T, the following condition should also hold: S's(Tq^^) < 
7x lOVVVd"^- If fs = 0.16 Hz and Ti > T (nonstationary), the following 
condition should also be true: 5,(Toi^) < 2 x lO^VVd"^- If /« = 0-32 Hz 
and Ti < T, the following condition should also hold: 6 x 10^ fiV'^f^^ < 
5,(Toi^) < 3 X lOVVVd"^- If fs = 0.32 Hz and Ti > T (nonstationary), 
the following condition should also be true: 3 x 10^/iV^/^^ < Ss{Tq^^) < 
3 X lOVVVrf-'. 

(4) Group IV. Necessary condition: /, < 0.32 Hz. If /, = 0.16 Hz, the fol- 
lowing condition should also hold: ^.(Tq-/) > 7xl0^fiV^f-\ If/, = 0.32 
Hz, the following condition should also hold: S's(Tq^^) > 3 x lO^yuV^/J'^ 

It can be seen from Table |2] that the difference between this partitioning 
of 84 subjects and medical diagnoses made by interviewing the subjects is 
20%, given that the subjects in groups I and II are mapped to the "healthy" 
group of the medical diagnosis and groups II and IV to the "susceptible to 
schizophrenia" group. 

It is known that the signals produced in vivo by a human organism con- 
tain both a collection of low- frequency ( "resonance" ) components and higher- 
frequency stochastic variability [17lllHlll9l|50l|5H|52l[53]. We believe that 
the random components are responsible for the adjustment of the organ- 
ism to constantly changing external conditions, which also manifests itself 
as the rearrangement of resonance components. Our results imply that this 
general feature of in vivo systems can be seen in the analysis of frequency- 
phase synchronization. Our study shows that the degree of disruption in 
frequency-phase synchronization can be considered as a diagnostic sign in 
identifying possible pathological changes in the higher neural activity of chil- 
dren/adolescents. To quantify frequency phase-synchronization and its dis- 
ruption, we introduced a parameter fs, specific frequency characterizing the 
synchronization depth, and FNS parameters for stochastic signal compo- 
nents, respectively. It can be assumed that parameter a, which varied within 
the same ranges for subjects in groups I-III and had slightly higher values 
for group IV, is a measure of the randomness responsible for the adjust- 
ment of the organism to changing external conditions. On the other hand, 
the spikiness factor ^^(Tq^^), which varied in the analyzed signals up to 4-5 
orders of magnitude and was highest for subjects in group IV, appears to 
be related not just to the stochasticity associated with jump- and spike-like 
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irregularities, but also to the disruption of frequency-phase synchronization. 

6. Concluding remarks 

Our analysis demonstrates a high specificity of EEG signals and the pro- 
cesses of their space-time synchronization, which may be attributed to an 
individual-specific set of natural frequencies and their interferential contri- 
butions rearranging in time. These rearrangement processes in the brain are 
most likely associated with the organism's response to the action of external 
factors. The specificity of EEG signals is also observed in the behavior of 
stochastic components, particularly, in the values of ENS parameter a, the 
standard deviation of the value of measured dynamic variable from the slowly 
varying resonant (regular) component, and in the flicker-noise dependences 
of stochastic components. This high level of specificity of biomedical sig- 
nals makes the diagnosis of functional disruptions in organism's subsystems 
complex. 

We show that in the case of schizophrenia one can substantially reduce 
the effect of specific features of biomedical signals by analyzing the quantita- 
tive measures of frequency-phase synchronization between the EEG signals 
recorded at different cortex areas. According to our working hypothesis, the 
occurrence of a psychiatric pathology leads to a significant deviation of the 
frequency-phase synchronization level from it is normal value. As the mini- 
mum level of disruptions under normal conditions may be individual-specific, 
the estimation of "individual norms" should become a topic of high concern. 
This topic is related to a more general problem of developing the princi- 
ples for personalized medicine, an emerging medical model emphasizing the 
customization of healthcare. 

We believe that the proposed working hypothesis can be applied to the 
diagnosis of other psychiatric disorders as well as neurodegenerative disorders 
(Parkinson's, Alzheimer's, and Hungtington's diseases; progressive supranu- 
clear palsy; and epilepsy). This implies that the diagnosis of any of the above 
disorders requires at least two digitized biomedical signals measured in dif- 
ferent cortex areas at the same time. In the presented study, we analyzed 
the EEG signals from scalp-mounted electrodes at locations F^ and F4 for 
84 children/adolescents and were able to classify them into 4 groups based 
on frequency-phase synchronization levels, which appear to be related (judg- 
ing from the comparison with medical diagnoses) to the subjects' risk levels 
for developing schizophrenia in adult age. For other psychiatric and neu- 
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rodegenerative disorders, appropriate cortex sites should be selected based 
on preliminary analysis or available literature. The cross- correlation analysis 
presented above can also be used to estimate the effectiveness of a thera- 
peutic method by comparing the values of frequency-phase synchronization 
levels before and after the therapy. 
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Appendix A. Parameterization algorithm in discrete form 

Consider a time series Vd{k). The subscript d here and below is used 
to denote the discrete form of expressions. Let Nf be the number of points 
corresponding to the selected averaging interval T, M be the number of 
points used in estimating the autocorrelation function. In this case, the 
parameterization procedure can be written as follows: 
1. Calculate the arithmetic mean for the signal: 



1 

f^v = j;vJ2^Ak)- (A.l) 

* fc=i 

2. Subtract the arithmetic mean from the series Vd{k): 

Vd{k) = Vd{k)-^Iv. (A.2) 

3. Calculate the autocorrelation function for the series V^: 

Nt-p _ 

Mp) = jr-^ E ^'^(^) VAk + p),P = 0..M. (A3) 

* ^ fc=i 



27 



Let the autocorrelation interval M be T/4 (higher values of M will result in 
the loss of statistical information in estimating the autocorrelation function). 
To go from discrete form to the continuous one, one can use the following 
expression: p = Nit/T. 

4. Calculate the discrete cosine transform of the autocorrelation function: 

SM = MO) + ^d (M) (-1)'' + 2 J2 Mp) cos (^) , (A.4) 

p=i 

where q = 0..M. For q = 1..M — 1, S'rf(g) should be multiplied by 2, which is 
the standard procedure for discrete Fourier transforms to take into account 
the spectral values in the second half of the frequency range. Here, relations 
q = 2ff^^M and Sd{q) = S{f) x fa describe the equivalence between the 
discrete and continuous forms of power spectrum estimate. 

5. Calculate Ssd{0) as the average value of the power spectrum for the points 
1 and 2 (point 0, which corresponds to the zero frequency, is not used in 
calculating Ssd(fi)): 



6. Interpolate |5'rf(g)| given by Eq. (A.4) using the expression: 



Ssdiq) = ^^'a^l ^ (A.6) 

by the method of nonlinear least-square fitting to determine the values of 
parameters n and Tqi^. The fitting is done on the basis of a double logarithmic 
scale, dividing the entire series into a set of equal intervals. We used the trust- 
region algorithm for nonlinear square fitting, which is built in MATLAB v. 7 
or higher [M] . 
7. Separate out the resonant component: 

Srd (g) = 5"^ (g) - Ssd {q),q = 0..M. (A.7) 



8. Calculate the autocorrelation function for the resonant component as 
the inverse discrete cosine transform of Srd{q)- When q = 1..M — 1, divide 
Srciil) by 2 to take into account the spectral values in the second half of the 
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frequency range. Then calculate the inverse cosine transform: 

1 



Ad{p) 



{Srd{0) + Srd{M){-iy} 

npq'- 



2M 



2M 

2 ^ Srd{q) COS 
q=l 



M 



9. Calculate the difference moment for the resonant component: 

<^S(P) = 2 [^.,(0) - i,,d{p)] ,P = 0..M. 

(2) (2) 

The continuous equivalent of $^d (p) i^ ^r (t). 

10. Calculate the difference moment for the experimental series: 



(A.8) 



(A.9) 



.(2) 



^Tip) 



Nt-p 



Nt-p 

E 

fc=l >- 



Vdik)-V,{k + p) 



11. Calculate the difference moment for the random component: 



.(2) 



.(2) 



(2), 



^:M = ^\'ip)-^iJip)- 



(A. 10) 



(A.ll) 



(1') (1^ 

The continuous equivalent of ^IJiiip) is $es (t). 

12. Determine the parameters a, Hi, Ti^ by fitting ^J^iip) in Eq. (A.ll) to 
the interpolation expression of the anomalous diffusion type[ 



$(2)(p) = 2^2 X [1 - V-\Hi) X T{Hi,p/T, 



' sd 



Id 



(A.12) 



where T{s,x) = /exp(— t)t* ^dt, T{s) = r(s,0), using the same least-square 

X 

fitting method as in step 6. 



13. Calculate SsdiT^ilj) by Eq. (A. 6). 

14. After the values of all six ENS parameters - a, Tqi^, Tid, Hi, n, SgdiT^^^ 
- are determined, calculate the dimensional values for Tqi^, Tid, SsdiTQ^^: 
To = Tow X At, Ti = Tid X At, S^iT^i^) = Ssd{T-^\) x At, where At = f^^ . 

15. Calculate the relative error e# in the interpolation of difference moment 
^f{p): 



eg, 



IVl 

E 
p=i 


*f (p) - -^S (p) - 


- $iJ (l>) 




p=l 





X 100%. 



(A.13) 
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Here, the error is determined as the ratio of the difference of areas between 
the experimental structure function and the total interpolation function to 
the area of the experimental structure function. The areas are calculated by 
numerical integration using the rectangle method because the original series 
have a rather large number of points. The parameterization is successful if 
e$ < 10% 



Appendix B. Cross-correlation function in discrete form 



In discrete form, Eq. (17) is written as 



qijdirir.ne) 



Nt-nT-U[ne] \ne\ 

E [Vid{k) - V^d{k + Ur)] [Vjd{k + ne) - Vjd{k + ne + n,)] 

k=U[-ng]\ng\+l 



jNt-riT-Ulng] \ng\ lNt-nr~U[ng\ \ng\+ne 

E [V^d{k)-Vid{k + n^)f J Z [V,dik)-V,dik + nr)f 

k=U[-ne]\ng\+l y k=U[-ne]\ng\+l+ng 

(B.l) 
1, X > 0; 



■) 



where Ur = [r/At\, ne = [%/AtJ, U [x] , ^ ^ ^^ 
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